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A nonperturbative proof is given that the chirally-invariant quark propagator contains both 
particle and hole singularities with different dispersion relations. Mesons made of a quark and a 
hole will produce dilepton pairs at masses characteristic of the plasma and with a distinct energy 
dependence. 
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I. INTRODUCTION 



The central purpose of the experimental program at RHIC is to produce a quark-gluon plasma. Despite 
a great deal of theoretical effort there is no clear experimental signature of the quark-gluon plasma that will 
distinguish it from the very complicated hadron phase. The Theory Workshop on Relativistic Heavy Ion 
Collisions is intended both to refine previous ideas and to initiate discussion of new possibilities. It is in the 
latter spirit that this paper proposes investigation of a new experimental signal not previously considered. 
The present analysis of this signal still leaves unanswered questions that require further investigation by 
other theorists using other techniques. 

In the QGP at temperature T, a typical state |$) will have almost all the quark and antiquark modes 
with momentum less than T occupied and considerably fewer modes occupied that have momentum larger 
than T. If a quark or antiquark is added to this background the result is | $) or d> | . The new possibility 
is that removal of a quark gives a state 6|$) and the norm of this quark vacancy state is not small if the 
momentum is less than T. Similarly an antiquark vacancy is d\$>). Chiral symmetry plays an essential role 
in this argument because the hole states are only important if &|$) has an energy higher than |$) and this 
is only possible when the constituent quark mass is smaller than T . Thus approximate chiral symmetry is 
essential. This has two consequences: (1) Only u and d quarks are candidates and (2) When chiral symmetry 
is badly broken, i.e. the hadron phase, there are no hole states. 

The proof of hole states in Sec II is entirely non-perturbative. The end of Sec II contains references to 
perturbativc results. Sec III focusses on experimental consequences of the hole states, the most important 
of which is the annihilation of a quark and hole into a virtual photon that subsequently converts to a lepton 
pair. If this quark and hole are bound into a meson, there will be a resonance in the dilepton invariant mass 
spectrum with new features. 



II. QUARK HOLES IN THE EXACT THERMAL PROPAGATOR 



A. Singularities of the exact propagator 

It is remarkable that certain simple properties of the exact thermal quark propagator guarantee that the 
existence of low momentum hole states (k < T) in the chirally symmetric plasma. The inverse of the exact 
time-ordered thermal propagator must be of the form 

[S(u, k)] _1 = l0 A{u, k) - 7 • kB{u, k) (2.1) 
because of chiral invariance. The propagator itself is 

S(uj ,k)= 1 -^-±±+ 1 -^±±±. (2.2) 
v ; 2 A — B 2 A + B y ' 

At k = the propagator cannot depend on the direction of k. Therefore B(oj,0) — and consequently 
S(u,0) — 7o/A(w,0). Now comes the one dynamical input, namely that at k = there is a pole in 
the propagator at the thermal mass of the quark. Lattice calculations of temporal quark propagators have 
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measured this thermal mass [JD . The thermal mass must have a negative imaginary part because of damping. 
By dimensional analysis, both the real and imaginary parts are proportional to temperature. Denote this 
complex mass by mj. Then at k = 0, A(u>, 0) = (uj — rriT)a(uj) where a{u>) does not have a pole at uj — tut- 
At non-zero momentum A has the str uct ure A(uj, k) = (uj — tot)o(w) + A\(uj, k) where A\ is defined to vanish 
at k = 0. The first denominator in (2.2) is A — B = (uj — rriT)a(uj) + Ai(uj, k) — B(uj, k). Since A\ and B 
both vanish at k = 0, they are both small when k is chosen sufficiently small. Therefore A — B will vanish 



at uj = E p , where E p is equal to tot plus a small, fc-dependent correction: 



En 



m T + [—A 1 (m T , k) + B(rriT, k)]/a(rriT) + 



(k small) 



(2.3) 



For small k the other denominator A + B vanishes at uj = Eh where Eh is equal to plus a different small 
k correction because B occurs with opposite sign: 



E h « m T + [— Ai(m.T, k) — B(m Tl fc)] /ci(tot) + 



(k small) 



(2.4) 



The first energy will be called the particle energy and the second the hole energy. The justification for 
these names will come in the next section. It is easy to see that the inclusion of a bare quark mass in the 
propagator spoils these arguments. 

Knowing that A— B vanishes at uj = E p and A + B vanishes at uj = Eh, one can show that PCT invariance 
requires A — B to also vanish at uj = —Eh and A + B to also vanish at uj = —E p . The exact propagator 
therefore has four singularities and may be written 



S{uJ,k) = i(7o -7 - k) 



Zp{uj,k) Zh{—uJ,k) 



uj — E„ 



E h 



+ + 7- k) 



Zp{—uj,k) Zh(uJ,k) 



E n 



uj - E h 



(2.5) 



where Im E p < and Im E h < 0. 



It is difficult to deduce much about the numerator functions Z p and Zh- 



They are very complicated functions of uj and contain all the branch cuts and any other singularities of the 
exact propagator. The ±uj in the arguments of these functions are determined by PCT. From the previous 
argument we know that at zero momentum the residues of the two kinds of poles are equal: Z p (mT,0) = 
Z h (rriT,0). 



B. Operator meaning of the hole singularities 



The e xist ence of a pole at a tot in the zero-momentum propagator led to the distinct poles at ±E p and 
±Eh in (2^). The next task is to figure out what these poles mean in terms of quark creation and destruction 
operators. At time t=0 the exact field operator for massless quarks in a finite volume V can be expanded as 



^(0, x) = — Y^{u a (k, s)b S s e lk - 3 + v a (k, s)dl 

k.s 



—ik-a 



(2.6) 



The field operator at any other time can be expressed as ip(t,x) = exp(iHt)ip(0, x) exp(— iHt). The exact 
time-ordered thermal propagator may be written in terms of two resolvents: 



S(uj,k) = -( 7o -7*.fc) Rp(uj,k)-Rh(-uJ,k) +^( l0 +l-k)(-Rp(-uj,k) + Rh{uJ,k)^ 



(2.7) 



R p (uj, k) = (bt 



1 



uj + ie — H 



4) R h {u,k) = {bK k 



H 



(2.8) 



where PCT was used to express antiquark resolvents in terms of quark resolvents. For compactness I have 
used the thermo field dynamics notatio n H = H — H and the brackets denote the expectation value in the 
TFD ground state Both resolvents ( |2.8| ) are analytic for Im uj > and have branch cuts along the real 
ui axis. 

The pole in (| 2.5|) a t uj = E p is in the lower half uj plane and has a characteristic spinor structure. 
Comparison with Q2.7[ ) shows that it must come from R p (uj,k) since Rh(-uj) is analytic for Im uj < 0. To 
understand what this means, recall that each state state |<I>) in the thermal average is an energy eigenstate. 
The existence of the pole at E p means roughly that the state &t|<I>) has an energy that is higher by E p . (More 
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precisely, b'\<&) is not an energy eigenstate but it has a non-zero overlap with an energy eigenstate that is 
higher by E p .) 

The new feature in (|2.5| ) is the pole at to = Eh in the lower-half to plane with the spinor structure 70 +7 ■ k. 
Comparison with ( |2.7|) shows that it can only come from Rh(u),k) since R p (—uj) is analytic for Im u> < 0. 
Therefore the state with one fewer particle than $ nevertheless has an energy that is higher by Eh- 
This state is called a quark hole or quark vacancy. 

The residue of the poles requires evaluating the numerator functions Z(oj,k) in ( |2.5| ) at the energy of 
the pole. At k — the holes have the same residue as the particles: Zh(rriT : 0) — Z p (rriT, 0). However at 
k > T the residue of the hole Zh(Eh, k) becomes very small. This means that high momentum hole states 
decouple from the field operator ip(x). It does not mean that high momentum hole states disappear, only 
that &|$) becomes a poor description. Since gluons couple to color currents -0A°7 A 'V' an d photons couple to 
the electromagnetic current ipj^ip, hole states with high momentum propagate have very small coupling to 
gluons and photons. 



C. Perturbative results on holes 



Nothing in the argument presented above is perturbative. However the observation that the thermal quark 
propagator contains two kinds of poles, one for quarks and one for holes, did arise in one-loop perturbative 
calculations (|-||] . The quark and hole states have been constructed to one- loop order in H . The dependence 
of the spectrum upon chiral symmetry has been studied in one-loop order J^j. The quark and hole states 
have been used in calculations of dilepton production [|| and of Higgs damping rates Jl(| . 



III. NEW MESONS USING HOLE STATES 



A. Binding of quarks and holes 

It is easy to envision how new mesons can arise in the chirally symmetric phase though it is difficult to 
prove their existence. If | $) is a color singlet state in the thermal ensemble then b k | $) is a color triplet with 
an energy higher by E p (k) after all the interactions between the and the plasma are included. Similarly 
6 / t'| ( I > ) is a color antitriplet with an energy higher by Eh{k') when all its interactions are included. A color 
singlet state b\pk > |$) should have an energy of approximately E p (k) + Eh(k') minus the color binding energy. 

In the chirally symmetric phase the color force is not confining but it is still attractive in color singlet 
channels. There are two different contributions to this binding: gluon exchange and instantons. The first 
may be thought of in terms of a Bcthe-Sa lpet er equation with a kernel for the exchanged gluons. Although 



the numerators of the quark propagator (2.5) have no mass term, it is the denominators that control the 
kinematics. This leads to differential operators that are approximately Klein-Gordon d 2 + m 2 -,, with m-r ~ 
150 MeV, acting on the coordinates of the Bethe-Salpeter wave function. This is analogous to the positronium 
problem and even non-relativistic analysis may be useful. 

A second source of color binding comes from non-trivial topologies, specifically the random instanton 
liquid advocated by Shuryak p"T[ . Here because of chiral symmetry the field operator satifies the equation 
"f^D^tp — with no chiral-symmetry breaking mass term. When is an instanton configuration the 
classical solutions of this produce a strong attraction in mesonic channels in which the fcrmions have the 
opposite chirality, e.g. attractive in the scalar channel but not in the vector channel. The lattice calculations 
of the Shuryak group confirm this picture. Note that although the thermal mass of the quarks is present, it 
is irrelevant for the instanton attraction. 



B. Quantum numbers of QQh mesons 

To determine the possible quantum numbers of the new mesons it is useful to examine the bilinear operators 
Jp(a;) = ip i (x)tf : jripj(x), where T is one of the sixteen Dirac matrices and t a are one of the four 2x2 matrices 
for isospin or 1 . The choice of T determines the value of J for the operator and its behavior under parity 
and charge conjugation. To determine the quark operator content of each Jp, we again use the expansion 
( p.6[ ) of the exact field operator. Integrating Jp(0, x) over three dimensional space then gives an operator 
describing a meson at rest. 
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For example, if T = 475 then the result is an operator for pseudoscalar 7r a and a mesons: 

J cfx^t^i^ix) = 2 l J2[bl(k,s)t^d](-k 1 s) - dii-lsWfr&s)]. (3.1) 

k,s 

This operator has the expected structure: tfcft for the creation of a quark-antiquark pair plus its PCT 
conjugate db. The pair has no orbital angular momentum, L=0, and no spin, S=0. Since the total angular 
momentum is J=0, it correctly describes the ir and r\ mesons as 1 Sq state of quark and antiquark. 
As an example of a meson operator that gives a quark paired with a hole, let T — 70. Then 

J d 3 ^(f)^. 7 o^(x) ^2^[bl(k, S )t^b J (k,s) + d. l (k, S )tt 1 d](k, S )]. (3.2) 

k,s 

If this operator acts on a typical state |$) containing many quarks and antiquarks it can create a new kind 
of meson at rest, b' (k)b(k)\<&) (or the PCT conjugate) with an energy approximately E p + Eh minus the 
binding energy. 

The chart below tabulates the quantum numbers of the possible mesons for all 16 Dirac matrices. The 
choice of T determines the value of J PC for that operator. For the last four operators in the chart, Jp 1 
contains both ordinary meson combinations tfcft and new combinations b^b. 
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The lightest QQ mesons are at the top of the chart; the heaviest at the bottom. Isotriplets are listed before 
isosinglets. The arrows in the chart indicate the opposite parity states that are degenerate in mass because 
of 811(2)^4 symmetry. The 811(2)^ multiplets have various dimensions: n + a is four dimensional; p + ai is 
six dimensional; w is a chiral singlet. The QQ mesons automatically have P = (— and C = {— 1) L+S . 
The QQh mesons in the second column have P = (— 1) L but C is not determined by the values of L and S. 



C. Decays and experimental signals of QQh mesons 

In the chirally invariant phase not all the QQ and QQh mesons will be bound. The combinations that are 
bound into isosinglets have a Zweig suppressed decay into two or more gluons depending upon J PC . Their 
decay into photons might be detectable. The isotriplets cannot decay into gluons but they could break up by 
collisions with free quarks in the plasma. Most importantly, the holes can only exist in the chirally invariant 
phase. Thus in a RHIC collision they can exist for a maximum of 5 or 10 fm/c, corresponding to a minimum 
width of 20 MeV to 40 MeV. 

In the random instanton liquid model p"l| ] the combinations that bind are those whose constituents have 
opposite chirality, viz T = 1, 175, 75707, 7o7- The last two of these would give QQh vector mesons with 
J PC = 1^ and 1 . The latter would produce a characteristic dilepton pair signal at its mass. It seems 
quite unlikely that this thermal meson mass would coincide with that of any known sources of dileptons. 
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The calculation of meson binding and of decay depends crucially on the residue function Zh{Eh, k). As 
mentioned earlier it is the same as the quark residue k = but at k > T it is negligible. This complica- 
tion helps distinguish the QQh dilepton signal from background. All dilepton rates fall exponentially like 
exp(— E/T) where E is the energy of the lepton pair. For normal sources such as p and u, the prefactor of 
this exponential is almost independent of E. For dileptons produced by QQh mesons the prefactor will fall 
rapidly with the dilepton energy E and be negligible at E <~ T. Thus in a plot of dilepton invariant mass 
that includes only low energy pairs such as E < 50 MeV or E < 100 MeV, the QQh signal is most likely to 
stand out. 
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